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ABSTRACT 

New images from the Hubble Space Telescope of the ERII radio galaxy Pictor A reveal a previously undis¬ 
covered tidal tail, as well as a number of jet knots coinciding with a known X-ray and radio jet. The tidal tail is 
approximately 5" wide (3 kpc projected), starting 18" (12 kpc) from the center of Pictor A, and extends more 
than 90" (60 kpc). The knots are part of a jet observed to be about 4' (160 kpc) long, extending to a bright 
hotspot. These images are the first optical detections of this jet, and by extracting knot flux densities through 
three filters we set constraints on emission models. While the radio and optical fiux densities are usually ex¬ 
plained by synchrotron emission, there are several emission mechanisms which might be used to explain the 
X-ray fiux densities. Our data rule out Doppler boosted inverse Compton scattering as a source of the high en¬ 
ergy emission. Instead, we find that the observed emission can be well described by synchrotron emission from 
electrons with a low energy index {p ^ 2) that dominates the radio band, while a high energy index (/? ^ 3) is 
needed for the X-ray band and the transition occurs in the optical/infrared band. This model is consistent with 
a continuous electron injection scenario. 

Subject headings: Galaxies: Active, Galaxies: Jets, X-Rays: Galaxies, Galaxies: individual (Pictor A) 


1. INTRODUCTION 

To date more than 40 extragalactic jets have detections in 
the IR-UV band although the majority of these are nearby, 
low-power PR I jets, where the speeds are slower, and pro¬ 
cesses and environments are different from those of quasars 
and ERII radio galaxies. While many PR II quasar jets have 
been detected in HST images, most are so distant that only a 
few point-like emission regions are found. Resolving struc¬ 
ture in the knots is rare; one example of resolved structure 
is 3C 273, which has an unusually bright jet in the opti- 
cal band and is close enough that HST resolves each knot 
(iMarshall et all l200ll: iJester et aT]|2006l) . The proximity of 
Pictor A and its PR II jet (located at a redshift of z = 0.035) 
makes it an appealing candidate for resolving small, dim fea¬ 
tures not detectable in more distant quasars. Observing these 
features would improve our understanding of the physical pro¬ 
cesses active within these jets. 

There is currently some uncertainty regarding the primary 
emission mechanisms active within these jets at kilopar- 
sec scales. Observations by the Hubble Space Telescope 
(HST) have shown that often the radio and X-ray spec¬ 
tra cannot be connected smoothly via a single synchrotron 
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model. This could be explained by a jet with relativistic 
bulk motion, where the high-energy emission originates from 
relativistically boosted inverse Compton scattering of cos- 
mic microwave background photons (IC/ CMB; proposed by 
iTavecchio et al.|[2005 ICelotti et al.ll2QQlh . Boosted IC/CMB 
emission has the attractive consequence of allowing radio ob¬ 
servations to help constrain the conditions and environment 
of the X-ray emitting particles; no additional electron popula¬ 
tions need to be invoked. When lacking better data, surveys 
have also assum ed boosted IC/CMB emission to constrain 
the je t geometry (IMarshall et al .1120051: iHarris & Krawczvnskil 
I20Q2L who tentatively classified Pictor A as a boosted jet). 
The IC/CMB model also was supported by data from a num- 
ber of early surveys. A majority of the systems studied by 
ISambruna et al.l (120041) were best explained by an IC/CMB 
model, and f ollow up observations of PKS 1136-135 and 
1150-f497 by ISambruna et al.l (120061) verified predictions of 
the IC/CMB model. 

Despite this early evidence, the data for a number of 
these jet s are no longer con sistent with a boosted IC/CMB 
model. iMever et alT(120151) used gamma ray observations 
to disprove the boosted IC/CMB model for PKS 0637- 
572, which had previous ly been the prototypic a l example of 
a boo st ed IC/CMB iet (ITavecchio et ^120001: ICelotti et"^ 
120011) . iMever & GeorganopoulosI (120141) also used gamma 
ray observations to rule out a relati vistically boosted j et in 
3C 273, building on the work of iJester et akl (|2007() and 
lUchivama et al.l (l2006b . The high energy emission f rom P KS 
1136-135, one of the jets used bv ISambruna et al.l (120061) to 
verify IC/CMB predictions, has since been found to have 
a high level of polarization whic h is incompatible with an 
IC/CMB model (ICara et al.l 120131) . The data of all of these 
systems were previously consistent with a boosted IC/CMB 
model (most notably PKS 0637-572), but better data have 
ruled out this explanation. Here we contribute evidence 
against such a model for Pictor A, which was also initially 
labeled a boosted IC/CMB jet. 

Pictor A is an PR II radio galaxy at a redshift of 0.035. (Por 













































2 


Gentry et al. 


Hq = 70.5 km/s/Mpc, at this redshift, 1" corresponds to about 
700 pc.) Pictor A’s pen cil-like X-ray jet wa s discovered in 
Chandra observations bv lWilson et al.l (l2001h . extending 1.'9 
from the core, oriented toward the resolved NW hotspot 4.'2 
from the core. The jet is barely visible in radio images due to 
its low brightness relative to the core and the extended lobe s 
(iPerlev et al.ll 19971: ISimkin et^ll999l: [Marshall et al.ll201Qh . 
Simkin et al. discovered Ha emission about 5" from the core 
along the direction of the jet to the NW hotspot but found 
no oth er optical emission from the jet. iHardcastle & Crosto^ 
(l2005h first reported a weakly-detected 1 keV X-ray counter¬ 
jet extending to the E hotspot; more recent X-ray data of M. J. 
H ardcastle et al . (2015 , in preparation) confirm this detection. 

IWilson et~^ (1200Ih . on the basis of a boosted IC/CMB 
model for the X-rays, showed that the magnetic field in the jet 
would need to be substantially below equipartition, ^2x10“^ 
G, assuming that the jet is Doppler boosted with a bulk Lorenz 
factor r = 2 - 6 at 6^ < 30° to the line of sight. However, 
IHardcastle & Crostonl (120051) argued that the well-constrained 
steep X-ray photon index favors a synchrotron model for the 
X-rays, in which case the magnetic field strength and boosting 
parameters are not well constrained. 

Examining Chandra X-ray images of the jet, we found 
evidence for fiares in the jet at 3 - 4.5a significance 
(iMarshall et al.ll20l'^ . The projected size of the jet is over 
150 kpc long (4') and about a kpc wide, so finding localized 
brightness changes on 1-year time scales was surprising. Us¬ 
ing the 5 GHz radio fiux density, the equipartition magnetic 
field is ^eq = 17/iG assuming no boosting. If the 1 keV X-rays 
are due to synchrotron emission, then the corresponding elec¬ 
trons have 7^7 x 10^ (for E = ^nieC^) and the synchrotron 
loss time is ~ 1200 yr, similar to the light-crossing time scale 
for a 1" sized source: 2000 yr, which is much longer than the 
observed variability time scale. One possible explanation for 
the variability is that the emission arises from a very small 
knot inside the jet, perhaps as small as 0."002, so HST obser¬ 
vations were needed to investigate whether there are features 
< 0."1 in size. In addition to investigating these variable re¬ 
gions, HST data offered the opportunity to isolate and com¬ 
pare SEDs of multiple components within the jet, which had 
never been done for Pictor A. 

Here, we report the first optical detection of Pictor A’s jet 
and its previously unknown tidal tail using the HST Wide 
Eield Camera 3 (yVFC3). These data were combined with 
data from the Chandra X-ray Observatory, as well as the Aus¬ 
tralia Telescope Compact Array, in order to obtain a multi¬ 
wavelength view of this jet. 

2. OBSERVATIONS AND DATA REDUCTION 

Chandra X-ray data (1 keV) were reduced and provided by 
M. J. Hardcastl e et al. (2015, in prep aration), with radio (5 
GHz) data from [Marshall et all ([201()[) . with the process fur¬ 
ther described in M. J. Hardcastle et al. (2015, in prepara¬ 
tion). In our published images we also make use of X-ray (1 
keV) c ontours produced from the data used bv [Marshall et al.[ 
dMinh . which are a subset of the observations by M. J. Hard¬ 
castle et al. (2015, in preparation). The rest of this section will 
focus on the optical data reduction. 

Three images with the WPC3 were obtained under the HST 
Guest Observer program (GO proposal ID 12261). The wide¬ 
band filters used were: P160W, P814W, and P475W (see Ta¬ 
ble [T]). The P160W image was taken with the infrared detec¬ 
tor of WPC3 (WFC3/IR), while the P814W and P475W im¬ 
ages were taken with the UV-visible detector (WFC3/UVIS). 


Table 1 

Hubble Observations of the Pictor A Jet 


Start Date 

Filter 

A" 

(nm) 

AA^ 

(nm) 

Exposure 

(s) 

2011 Feb 25 

F160W 

1536 

268 

2708 

2011 Feb 25 

F814W 

802 

154 

1200 

2011 Feb 25 

F475W 

477 

134 

1299 


^ Pivot wavele ngth, as dehned in the WFC3 handbook 
dDressellDOlSh . 

^ Passband rectangular w idth, as dehned in the WFC3 
handbook dDresselll2015l) . 

The raw images wer e reprocessed using DrizzlePac (see 
[Gonzaga et al l [201^ . providing images at 0."02 binning for 
UVIS images and 0."07 binning for the IR image. In the IR 
image, the host galaxy is quite bright and extensive, adding 
significant background to a large fraction of the image, which 
requires that it be removed carefully before searching for jet 
features. 

2.1. Galaxy Fitting 

The program gal fit ([Peng et al.[[201Q[) was used to gen¬ 
erate a model of the Pictor A host galaxy and also generate 
uncertainty maps. The generated models consisted of a PSP 
core (caused by the central AGN), along with 2 generalized 
Sersic profiles to account for the elliptical host galaxy. The 
PSP was extracted from a nearby bright star, visible to the 
north of Pictor A and marked in Pigure[T] The Sersic bulges 
are described by Sersic indices n and effective radii r^. We 
found that a positive Q parameter (“boxiness”) was neces¬ 
sary to describe the outer bulge, but the inner bulge could be 
sufficiently fit with Co = 0. 

In reality we expect the profile to be more complicated, as 
the tidal tail provides evidence of a recent merger. These 
galaxy fits were primarily used to remove the large-scale 
structure of the host galaxy, allow us to visually identify 
knots. We do not claim any knot detections within 15" of 
the core due the rapidly varying fit residuals and the crowd¬ 
ing of sources associated with the host galaxy. In general, we 
sought the simplest fit which would still allow for knot iden¬ 
tification and accurate fiux extraction. The galaxy fits do not 
provide precise photometry for Pictor A. Results can be seen 
in Tableland Pigures [HU The P814W and P475W images 
did not require galaxy subtraction for the knot regions to be 
clearly separated from the bulge light (see Pigure[3). PigureO 
features the jet, our primary region of interest. 

2.2. Knot identification 

Using our images (galaxy-subtracted in the case of the 
P160W image), jet knots were identified visually by looking 
for features coincident among radio, optical and X-ray bands. 
Pigure [3] shows the jet’s extent in all bands with X-ray and 
radio contours overlaid. 

Pour knots were identified, located 32", 43", 106" and 
112" from the core. We label these knots HST-32 through 
HST-112 based on their location, with filter sub-labels ap¬ 
pended as necessary (e.g. HST-32-P160W). These knots are 
marked on images Figures [TJO Figures [4ll5] show these knots 
in greater detail, with X-ray and radio contours overlaid. 

Overall, the IR image (F160W) is more sensitive than the 
UVIS images (F814W and F475W) due to the difference in 
detector sensitivities and exposure times. Consequently, many 
of the knots are most clearly detected in the F160W image. 
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Table 2 

Galaxy Fit Parameters 


Filter 

Sersic Component mcore^ 

(ADUpix-*) 

^bulge^ 


(") 

hjd 

(°) 

Co^ 

F160W 

outer bulge 17.5 

-3.0 

16.9 

2.3 

16.6 

0.6 

-91.3 

0.5 


inner bulge 


17.4 

1.6 

4.2 

0.8 

-94.6 

0 


^ Apparent (ST) magnitude of the PSF model of the core. 

^ Sky background, assumed to be uniform across the image. 

^ Apparent magnitude of the bulge component. 

^ Sersic index. 

^ Effective (half-light) radius. 

^ Ratio of the elliptical Sersic profile’s semi-minor and semi-major axes. 

§ Position angle of the elliptical Sersic profiles’ semi-major axis, defined to be zero if oriented north-south 
and increasing as the ellipse is rotated to the east. 

^ The “diskyness/boxiness” ratio, which is greater than zero for more box-like isophotes. 



RA (J2000) 


Figure 1. Image from the WFC3/IR, using filter F160W. Data were binned to a scale of O'.'Ol, and smoothed with a Gaussian with a FWHM of Of'l. A 
grayscale was applied with a logarithmic str etch function. Central core is marked with a green cross. Jet knots HST-32 through HST-112 are labeled, along with 
the locations of the X-ray flares reported bv IMarshall et al.l (1201011 and the star used for a PSF in the galaxy-fitting process. The tidal tail is marked by a tangent 
arrow from the innermost point at which it is detected. Bad pixel regions have been masked in white, resulting in the ovals within Pictor A. 
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Figure 2. Same as Figure[T] but now a galaxy-subtracted image from WFC3/IR, using filter F160W. After the subtraction, the jet knots and the innermost extent 
of the tidal tail are more readily visible. The subtraction resulted in strong systematic errors within 15" of the active nucleus. 
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We considered the possibility that what we observe are ac¬ 
tually unrelated background sources, rather than jet knots. 
FigureO shows a large number of optical background sources, 
which could be confused as jet emission. By looking for 
features coincident with observed X-ray flux within the gen- 
eral jet geometry, we reduce the probability of false positives. 
iMarshall et all (120101) previously ruled out the likelihood that 
unrelated X-ray background sources could provide flux den¬ 
sities similar to the ones we observe associated with these 
HST knots. While we are unable to conclusively determine 
whether this observed emission is produced by the jet, we 
conclude that the emission is likely from jet knots, given the 
correlation in locations and shapes between the X-ray, radio 
and optical features. 

The true extent of the optical component of the jet is prob¬ 
ably larger than the 4 knots mentioned in this work. Given 
the high density of background sources, we have tried to be 
conservative in what we claim is a knot. Other features are 
likely to become visible with deeper optical imaging. 

2.3. Measuring Jet Knots 

Gaussian 2D fits quantifled the location and size of each 
knot. Figure [6] gives a sample visualization of a Gaussian fit 
applied to HST-32-F160W, the best detected knot. Results for 
all knots and Alter images are given in TableO In cases where 
the fit was only marginally statistically signiflcant, we listed 
the fitted flux density, Sy, as an upper limit. Fit visualizations 
(similar to Figure for e ach knot and band can be found in 
the work of lGentr^ ( l20l4 Appendix A), but those visualiza¬ 
tions are not necessary for this analysis. 

These Gaussian fits then determined apertures, through 
which we could extract photometric flux densities. While we 
could have used the fitted flux density, using rectangular aper¬ 
tures allowed us to be consistent in how we extracted the ra¬ 
dio and X-ray flux densities. These apertures were defined 
to be 2" wide across the jet, and approximately 2 Gaussian 
fit FWHM long in the direction of the jet. (Background was 
subtracted using an adjacent region of identical size.) The lo¬ 
cation and size of each aperture can be found in Table [H In 
the case of the HST-112-F160W, the desired aperture extends 
beyond the edge of the image. Without knowing how many 
counts we are missing, we simply denote the counts observed 
as a lower limit. 

The Gaussian profiles also allowed us to estimate more ac¬ 
curate upper limits for the marginally detected flux densities. 
The shape of the F160W fit (clearly detected for all knots) 
was used to create artificial knots of various flux densities, 
which could be injected into the other Alter images. For each 
knot and each band, there was a critical flux density below 
which we could not reliably retrieve the artificial knot from 
the noise. We quote that critical value as our upper limits for 
the aperture flux densities. These are much more conservative 
upper limits than those suggested by the Gaussian fits, which 
only considered statistical error. We believe this injection and 
retrieval process more accurately reflects the systematic errors 
which dominate our uncertainties. 

This approach of injecting artificial 2D knots would not 
work for setting limits on the radio sources, as the radio data 
were already collapsed into a ID jet profile before we ex¬ 
tracted a flux density. To assess the significance of a radio 
detection, we estimated a background noise level of approxi¬ 
mately 1 mJy beam“^ (with a 5 arcsec^ beamsize) near the jet 
(which is significantly higher than the off-source noise level 


of 40 /iJy beam-1), HST-32 and HST-112 are clear 

detections at this noise level (the radio contours of Figures [3]- 
[5] begin at this noise floor), knots HST-43 and HST-106 are 
not clearly detected in the radio data, so we denote their flux 
densities as upper limits. This introduces uncertainties in the 
spectral models we will construct in Section [3] but these un¬ 
certainties do not affect our final conclusions. 

We also considered the possibility of foreground dust 
extinction affecting these flux densities, and have found 
that it plays a neg li gible role. Using the work of 
ISchlaflv & Finkbein^ (1201 Ih we predict no more than .15 
magnitudes (15%) of extinction in any of these optical bands 
(assuming a galactic column density). This is sufficiently 
small to neglect while interpreting our data. 

3. MODELING THE KNOT SPECTRA 

Using the aperture flux densities of each knot, we tested 
a number of non-thermal emission models. Three emis¬ 
sion mechanisms were plausible: synchrotron radiation from 
high energy electrons moving in a magnetic field; syn¬ 
chrotron self-Compton emission, where synchrotron photons 
are scattered off of high energy electrons (SSC); and in¬ 
verse Compton scattering of electrons that boosts photons 
from the cosmic microwave background to very high ener - 
gies (IC/CMB; iTavecchio et all 120001: Ic^otti et al J llomh . 
For more details on various models for the physical systems 
underlying these emission m echanisms, see the review of 
[Harris & Krawczvnskil (120061) . In all models we assume low 
frequency flux is provided by synchrotron emission; the label 
we attach to each model refers to the primary mode of X-ray 
emission, with radio and optical emission always produced by 
synchrotron emission. In all of the models discussed below, 
the SSC emission was significantly smaller than the IC/CMB 
emission, so we deem SSC emission negligible for the ob¬ 
served knots. The un-boosted IC/CMB emission was also too 
low to be shown on our plots, but we did analyze the amount 
of boosting that would be required for the IC/CMB emission 
to match the observed X-ray flux densities. 

In order to construct these models, we first assumed the ra¬ 
dio and optical flux densities to be from synchrotron radiation. 
We need the low energy spectral slope, but since we only have 
radio observations in one band, we need to use optical data 
to constrain the low energy spectral slope. From the radio- 
optical spectral slope, a (such that Sy oc z^"^), we can then 
determine the energy index, p = 2a + I, of the emitting elec¬ 
trons (assumed to have a power law distribution of energies 
Ne(E) oc E~P). Using those values for p (typically p ^ 2), the 
observed radio flux density allows us to determine the mag¬ 
netic field strength that would require t he minimal ene rgy in 
the jet, using the equations from IWorralll (120091) . The 
minimum energy system is approximately equivalent to a sys¬ 
tem in equipartition. Table [6] lists the results for p and 
for each knot. It should be noted that is only accurate to 
within a factor of 2 or so, gi ven the assump tions that could be 
made using the approach of lWorrallldM)^ . 

The radio-optical spectral indices differ significantly from 
the measured X-ray spectral indices - there must be a spec¬ 
tral break between the optical and X-ray frequencies, a soft¬ 
ening of the spectrum. The most likely option is that there 
is a break in the electron energy distribution itself. (For a 
continuous, but kinked distribution we will use the label bro¬ 
ken', for a distribution with a jump discontinuity we will use 
the label multiple populations.) A break in electron energies 
is predicted by the continuous injection model: if electrons 
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Figure 3. The overall F160W image (top) from Figure|2] cropped to feature the jet, along with identical regions from the F814W (middle) and F475W (bottom) 
images. The F160W image has been galaxy-subtracted, while the F814W and F475W images have not. All images have been smoothed identically, and have 
similar logarithmic stretch functions. X-ray contours (left, red) were produced using 7 levels, logarithmically spaced between 0.5-12.5 nJy arcsec“^. Radio 
contours (right, blue) were also produced using 7 logarithmically space d levels, ranging from 1.25-5 mJy beam“^. (These contour levels begin approximately at 
the noise level of each observation.) Knots are featured in Figures NT^ 


Table 3 

Gaussian Fit Results for Jet Knots 


Knot Label Bandpass 



/ 



PA' 


(^Jy) 

(") 

(") 

(") 

(") 

(°) 


HST-32 

F160W 

F814W 

F475W 

HST-43 

F160W 

F814W 

F475W 

HST-106 

F160W 

F814W 

F475W 

HST-112 

F160W 

F814W 

F475W 


2.2 ±0.9 

32.0 ±.5 

<0.18 


<0.01 


0.4 ±0.2 

42.7 ±.l 

<0.18 


<0.11 


0.2±0.1 

106.4 ±.8 

<0.10 


<0.01 


3.8±1.9 

112.1 ±.2 

2.2 ±.2 

111.9±.l 

1.2±.l 

112.0±.l 


0.1±.l 

3.3 ±.9 

0.6±.l 

0.5 ±.2 

0.4±.l 

0.4 ±.2 

0.4±.l 

1.2±.4 

0.36 ±.02 

0.85 ±.06 

0.39 ±.01 

1.19±.07 


0.9 ±.3 

1.1±.2 

0.3 ±.l 

-38.1 ±.9 

0.3 ±.2 

-52 ±3 

0.6 ±.2 

-.1±.3 

0.43 ± .03 

1.2±.l 

0.42 ±.03 

6.0±.l 


^ Position of Gaussian centroid relative to the core, measured along (x) or transverse to (y) a line with 
position angle of-79° (E of N) defining the jet. 

^ FWHMs of the 2D Gaussian fits, along (x) and transverse to (y) the jet. 

° Position angle of the 2D Gaussian major axis relative to the direction to the core, defined to be zero when 
aligned with the direction to the core and increasing when rotated counter-clockwise. 


are injected with a single power law of energies, the high¬ 
est energy electrons would experience the greatest radiative 
losses, resulting in a steady state distrib ution coinposed of a 
broken power law, with phigh ~ T^iow + 1 (lLongairll201lL Sec¬ 
tion 16.3.1). This is precisely our situation, if we assume all of 
the observed emission is coming from synchrotron processes. 
We infer pio^ ^ 2 for the lower-energy electrons responsible 
for the radio-optical emission, and phigh 3 for the higher- 


energy electrons responsible for the X-ray emission. 

This spectral break can also be used to gain information on 
the age of a knot. For an electron with a Lorentz factor 7, in a 
magnetic field of strength B we expect the electron to have a 
characteristic synchrotron cooling time scale of: 


fcooi = 2.4 • 10 ’ 



yr 


( 1 ) 
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4. Same as Figure^ but now cropped to feature knots HST-32 and HST-43. HST-32-F160W is shown in more detail in Figure[6] 








Figure 5. Same as Figured except now for knots HST-106 and HST-112. HST-112 lies on the edge of the IR detector and is cut off in the F160W band (top). 
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Figure 6. HST-32-F160W, along with a 2D Gaussian fit. The image has a 
color scale, inverted in intensity relative to the previous images. Red contours 
are from the observed optical data; green contours are from the 2D Gaussian 
fit to the data. Image has been rotated 11° clockwise, so the jet direction is 
horizontal, going left to right. A foreground star has been masked from the 
analysis (in white). 


Table 4 

Knot Aperture Regions 


Knot Label 

RA^ 

(hh:mm:ss) 

Dec^ 

(dd:mm:ss) 

(") 

Sy ' 
(") 

HST-32 

5:19:46.686 

-45:46:37.35 

6 

2 

HST-43 

5:19:45.711 

-45:46:34.82 

1.7 

2 

HST-106 

5:19:39.732 

-45:46:22.80 

2 

2 

HST-112 

5:19:39.197 

-45:46:21.79 

4 

2 


^ Box center (J2000) 

Box length, along the jet, which is defined to be in the 
direction -79°E of N 
° Box width, transverse to the jet 


(adapted from lLongaidl201 Ih . Evaluated at the break energy, 
Tbreak^ef*^, we get a characteristic time scale fbreak- For the 
continuous injection model, this characteristic time scale is 
expected to be equal to how long an injection event has been 
continuously occurring for a particular knot. For a model with 
multiple electron populations, this time scale corresponds to 
the maximum time that could have elapsed since the low en¬ 
ergy population was last accelerated. (The high energy pop¬ 
ulation must have been accelerated even more recently, given 
the shorter time scales at higher 7.) 

For each knot we constructed an SED consisting of syn¬ 
chrotron emission from a broken power law of electron ener¬ 
gies. The low energy indices were set by the radio and optical 
flux densities, while the high energy spectral indices were set 
by the observed X-ray spectral indices. If possible, a contin¬ 
uous but kinked electron energy distribution was used. If that 
was not possible (i.e. HST-112), we used 2 synchrotron spec¬ 
tra with differing normalization (a multiple population model, 
introducing a discontinuity in the spectrum). The results for 
each knot can be seen in Table [6| and Figures ITlfTOl 

These results can be contrasted with models which are 
used to explain the X-ray emission through IC/CMB emis¬ 
sion. Either through bulk Doppler boosting or lowered mag¬ 
netic field strengths, IC/CMB X-ray emission can be boosted 
relative to the lower-energy synchrotron emission. In the 
case of bulk Doppler boosting, we can use the X-ray fiux 
density predicted by an un-boosted model to infer the re¬ 
quired Doppler factor: 6}^^ = 5^,observed/5'^y,predicted, where 

6 = (r[l-f3 cos^])“^ for a bulk Foren z factor T and an angle 
0 between the jet and the line of sight (IWorrall & Birkinshawl 



I 2 OO 6 I) . This approach for X-ray emission only works for val¬ 
ues ofS < 10^; beyond that limit, assumptions in the scaling 
break down. All of our models are close to, or beyond this 
limit, indicating that unphysical levels of Doppler boosting 
would be required. A similar approach can be done for low¬ 
ered magnetic field strengths, following the scaling relation: 
(^req/^me) ^ ^observed/predicted, with the total energy 

(in the case of lowered magnetic field strengths) following 
the relation: F’req/F’me ~ (^req/^me)”^^’^^^- These scaling re¬ 
lations would require magnetic field strengths around 10^ be¬ 
low equipartition, with an energy roughly 10"^-^ times greater 
than the minimal energy. We consider such a configuration 
unphysical. 

Even if the amount of IC/CMB X-ray emission could be 
significantly boosted, we would expect that the spectral in¬ 
dex would be identical to that of the radio-optical band, 
whereas our data rule out this possibility. One might ar¬ 
gue that in addition to unphysical amounts of boosting re¬ 
quired, there might also be a break in the electron energy dis¬ 
tribution between the IC/CMB emitting electrons (lower en¬ 
ergy) and the synchrotron emitting electrons (higher energy), 
but that would imply a hardening of the energy distribution 
(Phigh ~ Plow “ !)• Such a configuration would be unusual— 
the opposite of the softening predicted by the continuous par- 
ti cle injection model—an d is difficult to explain. (The work 
of ISambruna et al.l (I2QQ4I) would have tentatively labeled Pic- 
tor As jet a synchrotron emitter on the grounds of our spectral 
data as well; see their Figure 4.) Between the unphysical lev¬ 
els of boosting and the unexplained hardening of the electron 
energy spectrum required for a boosted IC/CMB model, we 
rule out boosted IC/CMB as the primary emission mechanism 
for the X-ray component of Pictor As jet. 

4. TIDAL TAIL 

A tidal tail, clearly visible in Fig. [2l was an unanticipated 
find. This tail starts coincident with a radio source 18" (12 
kpc projected distance) north of Pictor A, sweeping to the 
west along a path 90" (60kpc) long, and ending 90" (60kpc) 
from the center of Pictor A. (Many galaxy subtraction variants 
were attempted, but none showed the tail closer than the radio 
source 18" from the center of Pictor A.) At its closest to Pictor 
A, the tidal tail is approximately 3" (2 kpc) wide, becoming 
at least 7" (5 kpc) wide at its farthest from Pictor A. The tail 
was only detected in our deepest image (the F160W data), but 
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Table 5 

Aperture Flux Density Results for Jet Knots 


Knot Label 

(Radio) 

(mJy) 

Su(P160W) 

(^Jy) 

Su(P^UW) 

(^Jy) 

Sui¥415W) 

(^Jy) 

^^(X-ray)^ 

(nJy) 

ax-ray"'’ 

HST-32 

3.2 

1.98 ±.05 

< 2.0 

<2.0 

1.68 ±.07 

0.96 ±.07 

HST-43 

<0.2 

.26 ±.01 

< 0.4 

<0.2 

0.14 ±.06 

. . . c 

HST-106 

<0.4 

.22 ±.01 

< 0.4 

<0.3 

0.39 ±.03 

1.12±.15 

HST-112 

1.2 

> 6.98^ 

2.47 ± .06 

1.56 ±.04 

0.27 ± .03 

0.74 ±.23 


^ From the data of M. J. Hardcastle et al. (2015, in preparation) 

^ Spectral index, a, such that Su oc . 

^ Insufficient counts to extract a spectral index; a = 2 assumed for flux density extraction. 
^ Aperture extended beyond the edge of the detector. 


Table 6 

Synchrotron Model Results^ 


Knot Label 

Plow 

Q^low 

Tbreak 

(10«) 

Phigh 

Cl^high 

Bme 

(liG) 

(ireq 

^break 

(yr) 

HST-32 

2.25 

0.62 

39 

2.98 

0.99 

33 

4-102 

576 

HST-43 

<2.25 

<0.62 

26 

3.15 

1.08 

49 

< 3 ■ 10^ 

392 

HST-106 

<2.45 

<0.72 

21 

3.24 

1.12 

88 

<2-W^ 

150 

HST-112 

1.90 

0.45 

21 

2.48 

0.74 

29 

3-10^ 

1388 


^ See §[3]for deflnitions of these parameters of the synchrotron model fits to each knot. 





Figure 10. Same as Figure [T] but for HST-112. A single broken power law 
was insufficient; a second power law (with different slope and normalization) 
was needed. 

we expect deeper observations would reveal it in other optical 
bands as well. 

The classification scheme of [Bridge et al.l (I201QI) suggests 
that this long tidal tail is well into the first passage of a merger. 
Dynamical estimates predict such a tail would take a few hun¬ 
dred million years to form, and would last another few hun¬ 
dred million years. This is certainly long enough for the for¬ 
mation of a jet at least 160 kpc long, suggesting a possible 
connection between Pictor A’s large galactic-scale jet and a 
recent merger with the coincident radio source. 

5. DISCUSSION & CONCLUSIONS 

Our data contain the first optical detection of Pictor A’s jet, 
but we only clearly detect at most a few knots rather than the 
entire length of the jet. Although these data are sparse, they 
allow us to construct physical models of the jet emission at 
multiple points along its extent. 

The first major result was these optical data give radio- 
optical spectral slopes, a, which differ significantly from the 
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X-ray spectral indices measured by M. J. Hardcastle et al. 
(2015, in preparation). This discrepancy rules out a single, 
unbroken synchrotron spectrum. It also makes the boosted 
IC/CMB hypothesis very unlikely. 

Our optical data show that two components of synchrotron 
emission are much more consistent with observations. The 
break between these components would results in a soften¬ 
ing of the electron energy distribution, around 7 ^ 10^. Such 
a softening would be a natural consequence of a continuous 
injection model, where electrons are constantly being accel¬ 
erated to an initial power law of energies, but radiative losses 
result in a steady state broken power law of energies. The in¬ 
jection time scales implied by these spectral breaks are quite 
short, in the hundreds of years, while the light-crossing time 
scale is at least 1000 yr for each knot. We caution that deeper 
optical observations could rule out this single population syn¬ 
chrotron model with a cooling break, at which point this spec¬ 
tral break time scale would no longer give the duration of an 
ongoing injection event. If the single population synchrotron 
model were ruled out, a multiple population model might be 
necessary, which was the case for knot HST-112. 

Ruling out a boosted IC/CMB model is be a key element 
to understanding the em i ssion o f Pictor A’s jet. The work 
of iHardcastle & CrostonI (l2005h provided independent evi¬ 
dence against an IC/CMB model for Pictor A’s jet emission, 
which has recently been confirmed by M. J. Hardcastle et al. 
(2015, in preparation). Our work, providing optical data, also 
reaches the conclusion that a boosted IC/CMB model for Pic¬ 
tor A’s jet emission does not match observed data. 

Still unconfir med is the poss i bility of unresolved regions 
hypothesized by [Marshall et al.l (1201 Ql) based on X-ray vari¬ 
ability data. Unfortunately, our images were insufficient for 
identifying unresolved regions definitively associated with 
knot emission. All currently detected optical knots are re¬ 
solved on scales similar to the width of the X-ray jet emis¬ 
sion. Deeper observatio ns might yet uncover the unresolved 
regions hypothesized bv lMarstiall et al.l (l2010h . 

Deeper observations would also help constrain our SEDs. 
While most of the UVIS knots are non-detections (HST-32 
through HST-106, in the F814W and F475W images), our de¬ 
tection of UVIS fiux from HST-112 allows us to put stronger 
constraints on emission models. Deeper observations would 
help us understand if the other SEDs are simply from a single 
population broken electron power law (the simplest explana¬ 
tion) or from two distinct electron populations (such as HST- 
112). If the underlying electron distributions are broken, but 
continuous (as expected for a single population continuous 
injection model), determining the break energy would help 
constrain the time scales of particle acceleration events along 
the length of the jet. If instead the electrons are from multi¬ 
ple populations (a broken and discontinuous distribution), the 
high energy population might also help us better understand 
how regions of Pictor A’s jet could flare on short time scales, 
while emission from the majority of the jet remained stable. 
Better data, and stronger limits, would help us answer these 
questions. 

Finally, we have noted that our data contain the first evi¬ 
dence for a tidal tail near Pictor A, but it has only been de¬ 
tected in the F160W band. Deeper observations in visible 
bands would allow the study of this tidal tail and its role in 
the dynamics of Pictor A. 
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